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Background: Olprinone decreases the cardiac preload and/or afterload because of its 
vasodilatory effect and increases myocardial contractility by inhibiting phosphodiesterase III. 
Purpose: The objective of this study was to characterize the population pharmacokinetics of 
olprinone after a single continuous infusion in healthy male volunteers. 
Methods: We used 500 plasma concentration data points collected from nine healthy male 
volunteers for the study. The population pharmacokinetic analysis was performed using the 
nonlinear mixed effect model (NONMEM®) software. 

Results: The time course of plasma concentration of olprinone was best described using 
a two-compartment model. The final pharmacokinetic parameters were total clearance 
(7.37 mL/minute/kg), distribution volume of the central compartment (134 mL/kg), intercompart- 
mental clearance (7.75 mL/minute/kg), and distribution volume of the peripheral compartment 
(275 mL/kg). The interindividual variability in the total clearance was 12.4%, and the residual 
error variability (exponential and additive) were 22.2% and 0.1 29 (standard deviation). The final 
pharmacokinetic model was assessed using a bootstrap method and visual predictive check. 
Conclusion: We developed a population pharmacokinetic model of olprinone in healthy male 
adults. The bootstrap method and visual predictive check showed that this model was appropriate. 
Our results might be used to develop the population pharmacokinetic model in patients. 
Keywords: phosphodiesterase III inhibitor, men, pharmacokinetic model 

Introduction 

Phosphodiesterase (PDE) III inhibitors induce myocardial contractility and peripheral 
vasodilation through accumulation of intracellular cyclic adenosine monophosphate 
(cAMP) in myocardial and vascular smooth muscle cells. 1 " 3 PDE III inhibitors increase 
cAMP levels without activating beta receptors and are thus useful for controlling 
acute exacerbation of chronic heart failure in patients resistant to catecholamines or 
those receiving beta blockers. In addition, PDE III inhibitors act without increasing 
myocardial oxygen consumption. 4 - 5 They are commonly used for perioperative circula- 
tion management, especially at the time of weaning from cardiopulmonary bypass after 
procedures such as aortocoronary bypass surgery and mitral valve replacement, which 
improves peripheral circulation and hemodynamic stability. 6 - 7 However, development of 
hypotension and arrhythmia after administration of such inhibitors should be carefully 
monitored. 1,7 Therefore, PDE III inhibitors are not established as a standard treatment 
for perioperative circulation management in the current clinical settings. Among PDE 
III inhibitors, olprinone exerts stronger vasodilatory effects, especially pulmonary 
vasodilatory effect and reduction in arterial afterload than milrinone. 3,8 ' 9 
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Target-controlled infusion (TCI) is an infusion system 
that allows physicians to adjust the infusion rate by using a 
pharmacokinetic model and to control drug concentration at 
the target site. 1011 Therefore, TCI of olprinone enables effec- 
tive and safe administration of olprinone. We performed a 
population pharmacokinetic analysis to develop a pharma- 
cokinetic model of olprinone by using plasma concentration 
data obtained from Phase I studies. 

Material and methods 

Subjects and studies 

In this study, we collected 530 plasma concentration data 
points of olprinone from nine healthy male volunteers 
(subject numbers 1, 2, 3, 4, 5, 6, 7, 8, and 9; total num- 
ber of subjects, 39) who participated in two clinical trials 
(Study I and II, Table 1). Missing data and concentrations 
below the lower limit of quantitation were excluded from 
the analysis. These two studies were performed to evaluate 
the safety, pharmacokinetics, and pharmacological effect 
after intravenous administration of olprinone. The eligible 
subjects were healthy adult male volunteers with no clinical 
abnormality in a medical interview, physical examination, 
electrocardiogram test, blood pressure measurement, blood 
test, and urine test. The study protocols were reviewed and 
approved by the Institutional Review Boards of each clini- 
cal study site. Written informed consent was obtained from 
each subject. 

Study design 

This clinical trial was conducted in two separate steps 
(Table 2). In the first step, seven single doses of olpri- 
none were administered as a 5-minute infusion (Study I). 



Table 1 Cumulative 


numbers of 


subjects and 


background 


factors 








Total 


Number of subjects 


9 






Number of measurements 


530 (500) a 








Study 1 


Study II 




Total number of subjects 


27 (26)' 


12(11)' 




Number of measurements 


312 (300) a 


218 (200)' 






Age 


Height 


Body 




(years) 


(cm) 


weight (kg) 


Mean 


27.2 


171.5 


64.1 


Standard deviation 


2.8 


4.4 


6.9 


Median 


26.0 


173.0 


61.5 


Minimum 


24.0 


163.2 


56.5 


Maximum 


32.0 


177.6 


75.0 



Notes: Subjects were excluded from this analysis. The numerical value in 
parentheses shows the numbers after exclusion of these data. 



Table 2 Study design 


Stage 


Dose (u.g/kg) 


Group 


Subjects 
(n) 


Subject 
numbers 


Study 1 










1 


5.0 


A 


4 


1,2, 3, and 4 


II 


10.0 


B 


4 


5, 6, 7, and 8 


III 


2.5 


B 


4 


5, 6, 7, and 8 


IV 


20.0 


A 


4 


1,2," 3, and 4 


V 


30.0 


B 


4 


5, 6, 7, and 8 


VI 


1.25 


C 


6 


1,5, and 9 


VII 


50.0 


A 


4 


1,2, 3, and 4 


Study II 










VIII 


61.25 (loading dose, 
20 u.g/kg; continuous 
infusion rate, 
0.25 |lg/kg/min) 


B 


4 


5, 6, 7, and 8 


IX 


122.50 (loading 
dose, 40 Hg/kg; 
continuous infusion 
rate, 0.5 ug/kg/min) 


A 


4 


1,2, 3, and 9 


X 


163.75 (loading 
dose, 40 Hg/kg; 
continuous infusion 
rate, 0.75 ug/kg/min) 




4 


1, 5,' 6, and 7 



Notes: Subjects in Stages I to VII received a 5-minute infusion. The group receiving 
1.25 flg/kg consisted of 6 subjects. The subjects were randomly divided into 
subgroups with three subjects in each subgroup. To compare the occurrence of 
hyperemia between the subgroups, three subjects were administered olprinone, 
and the other three were administered the placebo. Subjects in Stages VIII to X 
received a single 3-hour infusion. a Subject was excluded from the analysis; b B' group: 
one subject in the B group was replaced with a subject in the A group; c subject was 
excluded from the analysis. 

After Study I, three doses were administered as a 3 -hour 
infusion combination of a 165-minute maintenance infusion 
following a 1 5-minute loading dose (Study II). Each stage 
comprised four evaluation periods as follows: 1 ) a screening 
period; 2) a baseline period a 24-hour period before the start 
of the study treatment period during which baseline evalua- 
tions were made; 3) a study treatment period 24 hours after 
drug administration; and 4) a follow-up period 7 days after 
drug administration. The intervals between administrations 
were at least 7 days. 

Study I 

We randomized eight subjects to the two treatment groups 
(Group A and Group B). Subjects in Group A (subject num- 
bers 1, 2, 3, and 4) received 5, 20, and 50 |-ig/kg olprinone 
as a single continuous infusion over 5 minutes. Subjects in 
Group B (subject numbers 5, 6, 7, and 8) received 2.5, 10, 
and 30 jag/kg olprinone as a single continuous infusion over 
5 minutes. In addition, six subjects (Group C), including five 
who participated in Study I, were randomized in a double- 
blind manner to the either the 1.25 Ug/kg olprinone (subject 
numbers 1,5, and 9) or placebo subgroup. The treatments 
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were administered to each of the three subjects over 5 minutes. 
We collected 5 mL venous blood samples from each subject 
before the administration; at 5, 7.5, 10, 15, 30, and 45 minutes; 
and at 1 , 1 .5, 2, 3 , 4, 6, 8, and 24 hours after the administration 
to measure the plasma concentration of olprinone. 

Study II 

Subjects in Group B (subject numbers 5, 6, 7, and 8) received 
a loading dose at a rate of 1.33 ug/kg/minute for 15 minutes 
followed by continuous infusion at a rate of 0.25 |ig/kg/minute 
for 1 65 minutes; subjects in Group A (subject numbers 1 , 2, 3, 
and 9) received a loading dose at a rate of 2.67 |J.g/kg/minute 
for 15 minutes followed by continuous infusion at a 
0.5 |ig/kg/minute for 165 minutes; and subjects in Group B' 
(subject numbers 1, 5, 6, and 7) received a loading dose at a 
rate of 2.67 |ig/kg/minute for 1 5 minutes followed by continu- 
ous infusion at a rate of 0.75 |lg/kg/minute for 165 minutes. 
Group B' consists of three subjects from Group B and one 
subject from Group A of Study I. Group A consists of three 
subjects from Group A of Study I and a newly enrolled subject 
from Group C of Study I. We collected 5 mL venous blood 
samples from each subject before the administration, at 15 
and 30 minutes, and at 1, 2, 3, 3.04, 3.08, 3.16, 3.25, 3.5, 4, 
4.5, 5, 6, 7, 9, 11, and 24 hours after the administration to 
measure the plasma concentration of olprinone. 

Plasma concentration of olprinone 

The concentrations of unchanged olprinone in the plasma 
were measured using a sensitive and specific high- 
performance liquid chromatography (HPLC) method. To each 
1 mL of plasma, 1 00 \lL of internal standard solution (milri- 
none 10 ng/100 |lL methanol) and 1 mL of 0.2 M phosphate 
buffer solution (pH 5.5) were added and mixed on a vortex 
mixer. The mixture was applied on the top of a Bond Elut 
CI 8® column (Analytichem International, Harbor City, CA, 
USA). After washing the column with 2 mL H 2 0 and 2 mL 
H 2 0 containing 10% methanol, olprinone was eluted twice 
with 0.5 mL of methanol, and the effluent was evaporated to 
dryness at 50°C under a gentle stream of nitrogen gas. Then 
the dried residue was dissolved in 100 |iL of methanol and 
a 50 |lL portion was injected onto the HPLC column. The 
quantification limit was 0.1 ng/mL. 

This analytical method was validated for selectivity (no 
peak interfering with peak of olprinone), linearity (r=0.999), 
precision (intraday assay, 1.25% to 9.79%, n=4; interday 
assay, 3.24% at 4.90 ng/mL, n=3), recovery (89.2% to 125. 1% 
in 0.1-10 ng/mL, 92.5% to 107.3% in 10-100 ng/mL), and 
calibration curve (0.1-100 ng/mL). 



Development of a population 
pharmacokinetic model for olprinone 

The population pharmacokinetic analysis of olprinone 
was performed using the nonlinear mixed effect model 
(NONMEM®; ICON Development Solutions, Ellicott City, 
MD, USA) software package, version 7.2.0, with the PREDPP 
subroutines (ADVAN1, ADVAN3, or ADVAN1 1). The first- 
order conditional estimation with interaction method was 
used for estimation of the parameters. Bootstrap resampling 
was performed using the MULTTEST procedure of the SAS® 
9.2 software (SAS Institute Inc, Gary, NC, USA). 

We fitted one-, two-, or three-compartment models to the 
data when selecting the basic structural model. We exam- 
ined the interindividual variability for all pharmacokinetic 
parameters. The interindividual variability was calculated 
using an exponential error model. The residual variability was 
determined by using an exponential error model or a mixed 
error model. Dose and all pharmacokinetic parameters such 
as total clearance (CL), distribution volume of the central 
compartment (VI), intercompartmental clearance (Q or Q2 
and Q3) and distribution volume of the peripheral compart- 
ment ( V2 and V3 ) were adjusted for body weight. The model 
was selected by statistical evaluation of the objective function 
value difference (AOBJ) using % 2 -test, goodness of fit plots, 
and 95% confidence interval of parameter estimates. The 
95% confidence interval of parameter estimate was calculated 
using the point estimate +1.96 x standard error (SE), which 
was taken from the covariance step. Coefficient of variation 
(CV) or standard deviation (SD) of interindividual variability 
and residual error variability were calculated from square root 
of variance estimate. The significance level for the forward 
selection step and backward elimination step were 0.05 
(AOBJ <3.84) and 0.01 (AOBJ <6.63), respectively. 

Age as a covariate candidate was tested using the 
following formula: 

Pi = P + 0 Ba xBG (1) 

where Pi denotes the pharmacokinetic parameter estimates 
of subject i, P denotes the population estimates of the phar- 
macokinetic parameter, BG denotes the background factor 
of subjects, and G BG denotes the impact of background factor. 
The linear relation model was used because the range of the 
background factor was too narrow to use other models. 

Model validation 

The bootstrap method was used to evaluate the robustness and 
stability of the final model. 12 One thousand bootstrap data sets 
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were reconstructed by resampling the subjects from the original 
data set. The mean and SD of parameter estimates obtained 
from the bootstrap replications were compared to the final 
parameter estimates and SE obtained from the original data set. 
Successful run in NONMEM was defined as normal completion 
of both estimation and covariance steps. Moreover, a visual 
predictive check (VPC) was used to assess the final model. 

Results 

Subjects 

The demographic factors of the subject are summarized in 
Table 1. Nine subjects participated in this study. The total 
number of subjects in Study I was 27, and in Study II was 12. 
We collected 3 12 plasma concentration data points in Study I 
and 2 1 8 in Study II. Summary statistics (mean + SD) of age, 
height and body weight were 27.2+2.8 years, 171 .5+4.4 cm, 
and 64.1+6.9 kg, respectively. 

Plasma concentration-time profiles of olprinone for each 
dose (5, 10, and 2.5 |ig/kg in Study I) are shown in Figure 1. 
Plasma concentration of olprinone rapidly decreased after ter- 
mination of the continuous infusion. In subject numbers 2 and 5, 
plasma concentrations of olprinone were 0 at 1,440 minutes 



in the prior stage, and the dose intervals to the next stage were 
more than 7 days. Olprinone (0.053 ng/mL and 2.002 ng/mL) 
was detected in the plasma of subject numbers 2 and 5 before 
the administration. Therefore, these data (number 2 in Stage IV 
and number 5 in Stage X) were excluded from this analysis. 

Selection of the base model 

of population pharmacokinetic model 

The OBJ of one-, two-, and three-compartment models 
were 2,762.341, 1,835.064, and 1,268.120, respectively. 
The three-compartment model showed minimum OBJ. 
Generally, the OBJ decreases with an increase in the number 
of compartments. In the three-compartment model, the esti- 
mates of Q3 (0. 174 mL/minute/kg) and V3 (79.5 mL/kg) were 
smaller than the other estimates (CL, 7.34 mL/minute/kg and 
VI , 125 mL/kg; Q2, 8.74 mL/minute/kg and V2, 260 mL/kg), 
which indicated a minimum and clinically negligible effect 
of the third compartment on the time-concentration curve. 
Thus, we selected the two-compartment model, and the 
concentrations predicted using this model were mostly fitted 
to the observed data (Figure 2). Then, the residual error vari- 
ability was used for a mixed error model (exponential and 
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Figure I Plasma concentration-time profiles of different doses of olprinone (A— C). 

Notes: (A) O Subject I , • Subject 2, □ Subject 3, ■ Subject 4. (B) O Subject 5, • Subject 6, □ Subject 7, 



Subject 8. (C) O Subject 5, • Subject 6, □ Subject 7, ■ Subject 8. 
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A PREDvsDV B IPREDvsDV 




0 100 200 300 400 0 100 200 300 400 



PRED (ng/mL) IPRED (ng/mL) 

Figure 2 Goodness-of-fit plots (PRED versus DV, IPRED versus DV) for the final model. 

Notes: (A) Observed plasma concentration of olprinone versus predicted values. (B) Observed plasma concentration of olprinone versus individual predicted values. 
Abbreviations: DV, observed plasma concentration of olprinone; IPRED, individual predicted values; PRED, predicted values. 



additive model) (OBJ, 1,323.213). Therefore, the two- 
compartment model with the interindividual variability on 
CL was selected as a base model. 

Selection of the final population 
pharmacokinetic model 

Correlation coefficient between the post hoc estimate of phar- 
macokinetic parameter (CL), which was obtained from the 
base model, and the background factor of subjects (age) was 
0.240. In the forward selection step, age was not included as 
covariate on CL in the pharmacokinetic model (.P=0.363). The 
final model was therefore identical with the base model. 

Model evaluation 

Pharmacokinetic parameters from the final model are 
shown in Table 3. Estimates of CL, VI, Q, and V2 were 



Table 3 Population pharmacokinetic parameter estimates from 
the final model 





Mean ± SD 


%cv 


95% confidence 
interval 


CL (mL/min/kg) 


7.37±0.376 




6.63-8. 1 1 


VI (mL/kg) 


1 34+9.5 1 




1 15-153 


Q (mL/min/kg) 


7.7510.538 




6.70-8.80 


V2 (mL/kg) 


27519.97 




255-295 


co CL 2 


0.015310.00557 


12.4 


0.00438-0.0262 




0.049510.00548 


22.2 


0.0388-0.0602 


exponential 






^dtlve 


0.016610.00599 


0.129 (SD) 


0.00486-0.0283 



Abbreviations: to 2 , variance for interindividual variability; G 2 , variance for 
intraindividual variability; %CV, coefficient of variance; CL, total clearance; Q, 
intercompartmental clearance; SD, standard deviation; VI, distribution volume of 
the central compartment; V2, distribution volume of the peripheral compartment. 



7.37 mL/minute/kg, 134 mL/kg, 7.75 mL/minute/kg, and 
275 mL/kg, respectively. The CV of the interindividual vari- 
ability on the CL was 12.4%, and the residual error variability 
(exponential and additive) was 22.2% and 0.129 (SD). The 
relationships between the observed values and predicted val- 
ues (PRED) or individual predicted values of the final model 
are shown in Figure 2. The relationships between conditional 
weighted residuals 13 and time or PRED of the final model are 
shown in Figure 3 . The PRED obtained from the final model 
well described the observed data. 

Model validation 

The results of the bootstrap validation are shown in Table 4. 
The success rate was 99.3%, and the parameter estimates 
and SD of the final model were consistent with those of the 
bootstrap, which suggested good robustness and stability of 
the final model. Furthermore, the result of VPC is shown in 
Figure 4. The distribution of predictions simulated from the 
final model fitted well with those of the observations. 

Discussion 

We constructed a population pharmacokinetics model with 
full sampling data (average 13.6 points per subject) using 
nine healthy male volunteers. Tsubokawa et al 14 reported 
a population pharmacokinetics model for olprinone using 
26 patients undergoing cardiac surgery with cardiopul- 
monary bypass. In this study, the plasma concentration of 
olprinone was best described using a two-compartment 
model. The half-life of olprinone in the alpha and beta phases 
was 5.4 and 57.7 minutes, respectively. Tsubokawa et al 14 
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Time vs CWRES 




0 200 400 600 800 1,000 1,200 1,400 

Time (minutes) 

PRED vs CWRES 




100 150 200 

PRED (ng/mL) 

Figure 3 Goodness-of-fit plots (TIME versus CWRES, PRED versus CWRES) for 
the final model. 

Notes: (A) Conditional weighted residuals versus time. (B) Conditional weighted 
residuals versus the predicted values. 

Abbreviations: CWRES, conditional weighted residuals; PRED, predicted values. 



Table 4 Results of bootstrap analysis 


Successful 


993 






run* 


Final model 


Bootstrap 




Parameter 


Mean ± SE 


Mean ± SD 


95% confidence 








interval 


CL 


7.37±0.376 


7.39±0.376 


6.70-8. 1 7 


VI 


1 34±9.S I 


1 34±9.26 


I 1 7— 1 54 


Q 


7.7S±0.S38 


7.74±0.S52 


6.64-8.88 


V2 


27S±9.97 


275± 1 0.5 


25 1 -294 


W CL 2 


0.0IS3±0.00SS7 


0.0I39±0.00S39 


0.00486-0.0245 


exponential 


0.0495±0.00548 


0.049 1 ±0.00559 


0.0373-0.0592 




0.0I66±0.00S99 


0.0I69±0.006I8 


0.0077 1 -0.032 1 



Notes: a A successful run was defined as the number of normal completions of both 
the estimation step and covariance step of nonlinear mixed effect model among 
1,000 repeated analyses. Two runs failed at the estimation step. 
Abbreviations: &) 2 , variance for interindividual variability; a 2 , variance for 
intraindividual variability; CL, total clearance; Q, intercompartmental clearance; 
SD, standard deviation; SE, standard error; VI, distribution volume of the central 
compartment; V2, distribution volume of the peripheral compartment. 
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Figure 4 Visual predictive check in the final model for olprinone (Stage VIII, a 
loading dose at a rate of 1. 33 [Ig/kg/minute for 1 5 minutes followed by continuous 
infusion at a rate of 0.25 [Ig/kg/minute for I 65 minutes). 

Notes: (A) Linear plot; (B) semi-log plot. Open circles represent the data observed. 
Solid lines represent the 97.5th percentile and 2.5th percentile of the prediction 
interval. The dotted line represents the 50th percentile of the prediction interval. 



reported that CL, V 1 , and half-life were 7.13 mL/minute/kg, 
802 mL/kg, and 97.1 minutes, respectively. The half-life of 
olprinone in the beta phase in our study was shorter than that 
reported previously. Kimata et al 3 reported a positive correla- 
tion between the serum creatinine clearance and elimination 
half-life of olprinone. Tsubokawa et al 14 reported a mean 
creatinine clearance of 57 mL/minute, which was lower than 
the creatinine clearance in healthy male volunteers. Therefore, 
the difference in half-life was attributed to renal function. 

Clinical application of TCI infusion for drugs commonly 
used for anesthesia such as propofol, 1516 remifentanil, 17 and 
rocuronium 18 has been reported previously. However, to date, 
no study has described the clinical application of TCI infusion 
for PDE III inhibitors such as olprinone and milrinone. 

Olprinone is metabolized as a glucuronate conjugate in the 
liver. The plasma concentrations of olprinone that exert vaso- 
dilatory or cardiotonic effects in the extracted radial artery 
sections are 0.107 Limol/L (32.6 ng/mL) or 0.35 Limol/L 
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(106.7 ng/mL), respectively. 919 However, in the study by 
Arai et al 20 on off-pump coronary artery bypass, the blood 
concentration of olprinone 2 hours after administration at 
0.2 ug/kg/minute was 27 ng/mL. Subsequently, olprinone 
was administered at an infusion rate of 0. 1 ug/kg/minute, 
and the blood concentration was maintained over 1 5 ng/mL 
for 3 hours after returning to the ICU. Compared to the 
cardiac index after induction of anesthesia, that for 3 hours 
after returning to the ICU significantly increased. 20 Currently, 
PDE III inhibitors are not used as standard agents in peri- 
operative circulation management because of the possibility 
of development of adverse reactions such as hypotension 
and arrhythmia. After establishing the TCI infusion system, 
PDE III inhibitors may be used as standard agents in peri- 
operative circulation management. The TCI infusion system 
may be used in the same manner as the drugs such as propofol 
and remifentanil in the future; however, this system has not 
been validated for other drugs thus far. Further studies are 
required to develop a TCI system for olprinone. 

Our results showed that no background factors affected 
the pharmacokinetics of olprinone. However, our study 
included a limited number of subjects, limited study popu- 
lation, and limited range of covariates (age). Therefore, 
our results are not sufficient to conclude that no covariate 
affected the pharmacokinetics of olprinone. Further studies 
in patients are required to discuss the pharmacokinetics of 
olprinone. 

Conclusion 

We developed a population pharmacokinetic model of 
olprinone in healthy male volunteers. The final model was 
found to be appropriate using a bootstrap method and VPC. 
However, our study had limitations such as a small number 
of subjects, limited study population, and limited range of 
covariates. Thus, further studies in patients are required to 
determine the pharmacokinetics of olprinone. Our results 
might be useful to develop a population pharmacokinetic 
model in patients. 
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